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Abstract 

We discuss the discovery potential of the bound states of singlet quarks at LHC. 
We find that it is possible to discover bound states of singlet quarks at LHC with 
singlet quark masses up to 300 Gev for eg = | and up to 200 Gev for &q = — g. 
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The aim of this paper is the discussion of the discovery potential of the bound states 
of singlet quarks PJ|-[P] at LHC Singlet quarks are color-triplet fermions whose left 



and right chiral components are both singlets with respect to the SU(2) weak isospin 
gauge group. They can mix with the ordinary quarks and as a result of the mixing they 
decay into ordinary quarks. For small mixing e < O(0. 1) the lifetime of singlet quark is 
less than O(10) Mev so singlet quarks can form bound states. Note that as a result of 
mixing of singlet quarks with ordinary quarks tree level flavor-changing neutral currents 
are generated. Bounds from flavor changing neutral currents typically restrict the mixing 
parameter to be smaller than O(0. 1). So if singlet quarks exist they have to be rather 
longlived and have to form bound states. The phenomenology of bound states of new 



quarks has been discussed in refs. JL1J, [12|]. Recently it has been conjectured [13], [14| that 
singlet quarks can solve Rb - R c problem. 

In this paper we study a production and decays of _+ and 1 superheavy quarko- 
nium states. Other superheavy quarkonium states ++ , 1 ++ , 2 ++ have much smaller cross 
sections and they are practically unobservable [11, 12]. The main difference between our 
paper and refs. O, [T^] is that we use the LHC total energy i/i = 14 Tev and STEQ2L 



parton distributions. Besides refs. [|ll], [L2| consider the case of the 4 th generation dou- 
blet. Also we calculate the main backgrounds. The cross section for pseudoscalar 0^ + 
quarkonium production by two gluon fusion in pp collisions can be written in terms of 



gluonic decay width of + and the gluon distribution g(x, /x) in a proton as |TT|, [15 



a(pp ^ gg ^0 + 
7r 2 rr(0-+ -> gg) n 



8M 3 _ 



g(x,fi)g(x r, n)x l dx (12) 



M 2 _, 

with r = ° s . In nonrelativistic model the decay width of pseudoscalar quarkonium 
into two gluons is given by the formula |], 

r(o-+- 9S ) = ?f!^W)| 2 (i) 

Here Rs(0) denotes the radial wave function of the S state at r = 0. For the Coulomb 
potential we have |-Rs(0)| 2 = 4(|a s (Mg_ + )M -+) 3 . For other phenomenological potentials 



(Cornell potential, Richardson potential, Wisconsin potential) the value of the radial wave 
function at zero radius is higher [1], [Tl| so the use of the Coulomb potential gives the lowest 
cross section of _+ pseudoscalar quarkonium production. In this paper we shall use the 
Coulomb potential. The best way to look for the pseudoscalar quarkonium is through 
its decay into two photons. The reaction pp — > 0~ + — > 77 is very similar to the famous 
reaction pp — > Higgs — > 77 which is supposed to be the main reaction for the Higgs boson 
discovery at LHC. The decay width of pseudoscalar quarkonium into two photons is 

r(o- + -77) = ^^l^(o)| 2 (2) 

Here cq is the charge of the singlet quark Q in units of the proton charge. We have 
calculated the cross sections for the pseudoscalar quarkonium production using the parton 



distributions STEQ2L of ref.[16]. In our calculations we have used PYTHIA 5.7 and 



JETSET 7.4 generators ||17|| . We took the value of the renormalization point fi equal 



to the mass M -+ of the pseudoscalar quarkonium. We put strong coupling constant 
a s (M|) = 0.120, effective electromagnetic coupling constant a = a(M|) = and top 
quark mass m% = 175 Gev. We have checked also that the variation of the renormalization 
point fi in the interval 0.5M -+ — 2M -+ leads to the variation of the cross sections less 
than 30 percent. In our estimates we assumed the integral luminosity / L = 10 5 (p6) _1 
and the total energy y/s = lATev. The main background for the reaction pp — ^ h — > 77 
comes from: 

1. prompt diphoton production from quark annihilation and gluon fusion diagrams 
and bremsstrahlung from the outgoing quark line in the QCD Compton diagrams 

2. background from jets, where an electromagnetic energy deposit originates from the 
decay of neutral hadrons in a jet or from 1 jet + 1 prompt photon. 

The jet background is reduced |1| by imposing an isolation cut, which also reduces the 



bremsstrahlung background. The photon is defined an isolated [|18] if there is no charged 
track or electromagnetic shower with a transverse momentum greater than 2.5 Gev within 
a region AR > 0.3 around it. It is assumed that the jet background is reduced to an 
insignificant level (< 10 percent) by the combination of isolation and 7r° rejection cuts 



. For CMS detector [18 an efficiency of 64 percent was assumed for reconstruction of 



each photon (i.e. 41 percent per event). We assumed that the accuracy of the restoration 
of the diphoton invariant mass is 1 percent that is conservative estimate for CMS detector 
where the photon energy resolution is assumed to be AE/E = 0.02/ E ' 5 © 0.005 © 0.2/ E 
in the barrel and AE/E = 0.05/E ' 5 © 0.005 © 0.2/E in the endcap, where there is a 
preshower detector. In our calculations we have used the following cuts for photons: 

\v\ < 2-5, \p t \ > 25Gev 

The results of our calculations are presented in table 1. Our main conclusion is that at 
LHC for integral luminosity L = 10 5 pb^ 1 it would be possible to discover the pseudoscalar 
quarkonium at the 5a level for e q = — | with singlet quark masses uiq up to 100 Gev and 
for eQ — \ with singlet quark masses up to 300 Gev. 

For the estimate of the vector quarkonium 1 production we use the cross section 



for the subprocess gg — > 1 g equal to |JTT], [12 



where 



97T 2 S 

*(99 - 1-9) = 8MM7r2 _ 9) r (l- ^ 999)1^) , 0) 



w . 2 ,x+l 2xln(x) . 2{x-l) Unix) 

= -( — - - - — ^ + ; — f 2 + ^ 3 ( 4 ) 

x z x — 1 (x — 1) x(x + l) z [X + \) 6 



The decay width of vector quarkonium into 3 gluons is determined by the formula ]1|, [12 



40(tt 2 - 9)a 3 s 
81tt Mf. 



999)= \_„3 ; s \Rs(0)\ 2 (5) 



The most promising decay mode for the detection of vector quarkonium is 1 — > fi + fi~ 
or 1 — > e + e _ . The main background for the process pp — > 1 + ... — > [i + + ... is the 
Drell-Yan process. We have calculated the corresponding cross sections for singlet quarks 
with electric charge Eq — — | (analog of (i^-quark) and for singlet quark with electric 
charge eQ = | (analog of li^-quark). The results of our calculations are presented in table 
2. Our main conclusion is that for bq — — | it would be possible to discover singlet quarks 
with the mass up to 200 Gev and for eQ = | it would be possible to discover singlet quarks 
with the mass up to 250 Gev. 



It should be noted that an account of the production of the radial excitations of 0~ + 
and 1 quarkoniums leads to the increase of the signal cross section approximately by 
20 percent [[II] . 

To conclude, in this note we have studied the perspectives of the discovery of the bound 
states of singlet quarks, namely pseudoscalar quarkonium 0~ + and vector quarkonium 1 
at LHC The best way to detect such bound states is the measurement of the diphoton 
and dilepton invariant masses. LHC will be able to discover the superheavy quarkoniums 
with singlet quark masses up to 300 Gev for eg = | and and up to 200 Gev for tQ = — ~. 

I am indebted to the collaborators of the INR theoretical department for discussions 
and critical comments. The research described in this publication was made possible in 
part by JSPS program on Japan-FSU Scientists Collaboration. 
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Table 1. The cross sections and branchings for the process pp — > 0~ + + ... — > 77 + .... 
Here: 

1. (7 = er(pp — > 0~ + — > (?(?) • r _1 (0 _+ — > (7(7) in pb ■ (Mev)^ 1 . 

2. a cut = a(pp -> 0-+ -> ^||77 S | > 2.5, |p tr , ff | > 25Get;) •r~ 1 (0-+ -> <?<?) in pb- (Mev)- 1 

3. T = r(0 _+ -> in Mev. 

4. cr 7jCMt = cr(pp -> (T+ -> 77||?7 7 | < 2.5, |p tr;7 | > 25Get;) in /& • lO^ 1 for e Q = -|. 

5. o-Bacfc = ^Bacfcd^l < 2.5, |p tr , 7 | > 25Gev) / dm^ ■ 10~ 2 M -+ in /&. 

6 ' Ll = V^Bacfc ' ^ 5 = a ^fiut-L, N Back = a Back- L, L = 10 5 -pfc _1 , fc = 0.64, e Q = -§. 

7. L 2 = ^.g^, N s = 16k ■ <7 7 , cut ■ L, N BacK = a Back ■ L, L = 10 5 • pb~\ k = 0.64, 

e Q = l 

8. Mq is the mass of the singlet quark Q. 



M Q (Gev) 


a 


&cut 


r 




O'Back 


L x 


L 2 


100 


22.4 


16.6 


1.8 


75 


64 


6.0 


96 


125 


7.3 


5.6 


2.0 


29 


35 


3.2 


50 


150 


3.9 


3.1 


2.1 


18 


21 


2.6 


41 


175 


1.9 


1.6 


2.3 


10.7 


12 


2.0 


31 


200 


1.0 


0.83 


2.3 


5.7 


8.4 


1.3 


20 


225 


0.40 


0.34 


2.4 


2.7 


6.3 


0.69 


12 


250 


0.31 


0.27 


2.5 


2.2 


4.4 


0.67 


11 


275 


0.19 


0.17 


2.6 


1.4 


3.4 


0.50 


8.0 


300 


0.12 


0.11 


2.7 


0.98 


2.5 


0.39 


6.3 


350 


0.053 


0.047 


2.8 


0.46 


1.5 


0.25 


3.9 


400 


0.026 


0.023 


3 


0.25 


0.96 


0.16 


2.6 
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Table 2. The cross sections and branchings for the process pp — ■> 1 +... — > + 
Here: 

1- cr cu t,i = -> 1 + ••• + -IM ^ 2 - 5 ' |Ptr, M l > 5Get;) in fb for Q = -|. 

2. o- cuti2 = cr(pp -> 1" + ... -> + ...||^| < 2.5, |p trj/1 | > 5Gew in fb for e Q = §. 

3. VBack = da DY (pp -> + ...||^| < 2.5, |p tr J > hGev)/dm^ ■ M x — • (2 • 10 -2 ). 

5 - ^ 2 = T^T' ^ = °" cut > 2 ' L ' JVb «* = ° Back ' L ' L = !0 5 P& _1 , e Q = |. 
6. Mq is the mass of singlet quark Q. 



Mg(Gev) 




&cut,2 


O'Back 


u 


L 2 


100 


48 


80 


176 


36 


60 


125 


17 


29 


78 


19 


33 


150 


10 


17 


44 


15 


26 


175 


5.1 


8.5 


26 


10 


17 


200 


2.8 


4.6 


17 


6.8 


11 


225 


1.1 


1.9 


12 


3.1 


5.4 


250 


0.91 


1.5 


8 


2.8 


5.1 
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